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I. Introduction 

Microgravity processes must rely on mechanisms other than bouyancy to move bubbles or droplets 
from one region to another in a continuous liquid phase. One suggested method is thermocapillary 
migration in which a temperature gradient is applied to the continuous phase 1 . When a fluid particle 
contacts this gradient, one pole of the particle becomes warmer than the opposing pole. The intcrfacial 
tension between the drop or bubble phase and the continuous phase usually decreases with temperature. 
Thus the cooler pole is of higher interfacial tension than the warmer pole, and the interface is tugged in the 
direction of the cooler end. This thcrmocapiliary or thermally induced Marangoni surface stress causes a 
fluid streaming in the continuous phase from which develops a viscous shear traction and pressure gradient 
which together propel the particle in the direction of the w armer fluid. 

A significant and as yet unresolved impediment to the use of thermocapillary migration to direct 
bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid 
particle surface of surface active impurities dissolved in the continuous or (if the particle is a liquid) droplet 
phases. Surfactant adsorbs onto the surface of a moving fluid particle, where it is convected by the surface 
flow to the particle's trailing end Accumulation at the back end causes kinetic desorption into the bulk 
sublayer, and the sublayer concentration increases above the value far from the interface. This difference 
gives rise to a diffusive flux away from the trailing end. Similarly at the front end kinetic adsorption 
occurs from the sublayer since the front surface is swept clean of surfactant. The sublayer concentration 
adjacent to the leading end of the particle decreases creating a bulk diffusive flux from the bulk to the front 
end. Eventually a steady state develops: In this state, the surface concentration at the back end has 
increased to the point where the desorption rate, proportional to the difference between the surface and 
sublayer concentration, balances the convective rate. In addition, the sublayer concentration has increased 
sufficiently so that the diffusive flux away from the particle surface, proportional to the difference between 
the sublayer and far field concentration, balances the kinetic desorption. At the front end, the surface 
concentration becomes reduced enough so that kinetic adsorption balances convection, and the sublayer 
concentration becomes reduced enough so that diffusion to the surface balances adsorption Consequently, 
in this steady state the surface concentration is considerably higher at the rear than at the front of the 
particle, and the intcrfacial tension is lower at the back relative to the front. This intcrfacial tension 
difference creates a surfactant Marangoni stress along the surface as the front end tugs at the rear. The 
direction of this surface stress is opposite to that of the surface flow caused by the thermocapillary driven 
motion of the particle, and thus the adsorption of surfactant onto the particle interface acts to reduce the 
surface flow and hinder the intcrfacial mobility. The less mobile an interface, the more drag is exerted by 
the continuous phase on the particle as it moves through the medium, and the smaller is the 
thcrmocapiliary migration velocity. 

This reduction in surface mobility and migration velocity due to retarding Marangoni gradients 
caused by the convective partitioning of surfactant has been studied extensively for the buoyancy driven 



motion in which fluid inertial is negligible (low Reynolds number) and the fluid particles take the shape of 
spheres (small Weber number). For this case, the terminal velocity becomes reduced from the clean, 
mobile surface value (the Hadamard-Rybczynski velocity, U\\\t= pga 2 /3p (for a bubble) where p and p are 
the continuous phase density and viscosity, respectively, g the gravitational acceleration, and a the particle 
radius) to the Stokes velocity (Ust= 2pga 2 /9p) of a sphere when the interface mobility is completely 
arrested. When either the kinetic or diffusive transport is slow relative to interfacial convection, surfactant 
collects at the trailing pole in a stagnant cap of surfactant, and the terminal velocity is a function of the cap 
angle which is determined by the steady amount of surfactant adsorbed. 2 ? Theoretical studies for finite 
kinetic and diffusive transport have also been studied 61 °. The effect of surfactant adsorption on 
thcrmocapillary motion has only recently been studied: Kim and Subramanian 1 1 12 and Nadim and Borhan 1 " 
have theoretically examined the reduction in thermocapillary migration velocity due to surfactant 
adsorption, and find that, for the same amount of surfactant adsorbed onto the surface, the reduction 
(relative to the v<*|v$fbr a completely mobile interface) is much greater for thermocapillary driven movement 
than it is for buoyancy driven motion. The reason, as they point out, is that in thermocapillary migration 
the strength of the streaming flow around the particle which is driven by the surface tension gradient, and 
which is responsible for propelling the particle forward, is directly proportional to the interface mobility . 
The less mobile the surface, the more reduced is the streaming flow due to a fixed surface tension 
gradient, and the smaller is the particle velocity . Thus thcrmocapillary migration is very sensitive to the 
surfactant adsorption, and this adsorption, even arising from surfactants dissolved in trace amounts, can 
reduce the thermocapillary movement to near zero. 

Most experimental evidence in tests both on earth and under microgravity conditions indicate that 
it is difficult to achieve significant thcrmocapillary' migrations in agreement with the Young ct al value, 
and many studies have attributed the reduced migrations to the retarding effect of surfactant impurities 
unavoidably present in the bulk phases (sec the review article of Subramanian 14 ). In addition, Barton and 
Subramanian 1 ^ demonstrated directly the retardation by the intentional addition of surfactant to a liquid 
phase in which droplets were moving by thermocapillarity. To date the problem of retardation of 
thcrmocapillary driven particle motion due to surfactant impurities remains unresolved. Thcrmocapillary 
motions in agreement with the expression of Young et al 1 can only be achieved in systems in which 
extreme precautions have been taken to remove impurities (as, for example, in Barton and Subramanian’ s' s 
experiments using ethyl salicylate drops in ethylene glycol). 

In this paper, we provide a theoretical basis for remobilizing surfactant retarded fluid particle 
interfaces in an effort to make more viable the use of thcrmocapillary' migrations for the management of 
bubbles and drops in microgravity The retarding Marangoni stresses arise because the rate of convection 
of surfactant to the trailing pole is much larger than either the kinetic or diffusive flux away from the 
particle surface. The scale for the convective flux is T 0 Ua, where U is the terminal velocity , a (as before) 
the radius and T 0 is the surface concentration in equilibrium with the bulk concentration C G far from the 
particle. In this study we will use Langmuir relations to describe surfactant exchange; thus the adsorption 
rate is pCstT^-T) and the desorption rate is aT, C s the surfactant concentration adjacent to the interface 
(the sublayer concentration), V f is the maximum packing density , and a and P arc kinetic coefficients At 
equilibrium, the sublay er concentration is equal to C 0 and the equilibrium surface concentration To is 

T CJb 
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r„ i +cjb (1) 

where b is an adsorption parameter equal toc^. The scale for the diffusive flux is [DCo/aJ a 2 , where D is 


the bulk diffusion coefficient We assume surfactant concentrations arc below the critical micelle 


concentration, so no surfactant aggregates arc present in the bulk liquid. The scale for the kinetic 
desorptive flux isoCTo a 2 Retarding surfactant gradients do not develop when 
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where x = ba IV tr and Pe is the Pcclet number Pc=Ua/D, because when these inequalities arc met, kinetic and 
diffusive transport arc sufficient to maintain the surface concentration uniform despite the action of 
convective redistribution. The diffusion criteria in (2) illustrates the reason impurity or small 
concentrations of surfactant tend to retard interfaces: The diffusion coefficient of surfactant monomers is of 
the order of 1 0 6 cm 2 /sec, and b / T^ is of the order .01-1, so for typical values of a and U of 10 3 cm<a< 1 0** 
cm and 10 3 cm/scc < U< 1 cm/scc, y/Pe is not larger than I. Thus diffusion alone reduces the surface 
mobility. However, the criteria in (2) suggest a method for remobilizing interfaces which have been 
retarded by the adsorption of a surfactant impurity. Select a surfactant (which we term a remobilizing 
surfactant) whose kinetic desorption rate constant is larger than the convective rate U/a. There arc not 
many reported measurements of surfactant kinetic rate constants, but the small amount of data derived 
from dynamic surface tension measurements indicates that a can be as large as 10 2 see' 1 , so that for 
lJ/a< 1 0“ see' 1 , surfactants whose desorptive rate is fast enough so that (2a) is satisfied can be identified 
The second criteria can be satisfied if the bulk concentration is large enough; typically (,/b can be as large 
as 1 O' or more, before Micellar aggregation initiates. When the bulk concentration of the remobilizing 
surfactant is high, the impurity docs not have an opportunity to adsorb onto the surface because of the 
faster transport of the remobilizing surfactant Thus the interface consists only of the remobilizing 
surfactant, which keeps the interface mobile by maintaining a uniform surface concentration. 

To verily the above scaling arguments, we undertake numerical simulations of the effect of 
increasing the bulk concentration on the steady velocity of a gas bubble in a continuos liquid phase. In this 
first effort, we assume that buoyancy drives the bubble motion, and that surface tension forces arc larger 
than inertial and viscous forces (small Weber and capillary numbers) so that the bubble retains its spherical 
shape as it moves. The kinetic exchange is assumed to be fast enough to insure that the surface and the 
sublayer arc in equilibrium, as we have observed that surfactants can have fast kinetic regimes, and Pcclet 
numbers for the bulk surfactant mass transfer arc assumed of order one, as observed above. In addition, we 
retain the effects of fluid inertia by developing order one Reynolds number solutions since the Reynolds 
number is not small for bubble motions in the most common of continuous phases, w ater. In the follow ing 
sections we detail the Formulation and the Numerical Solution Algorithim (See. II ), and the Results (See. 
HI). 


II. Formulation and Solution Algorithim 

We consider the axisymmctric, steady motion of a spherical bubble in an unbounded, Newtonian 
incompressible liquid in a spherical coordinate sy stem (r,0) with the bubble fixed, and the coordinate angle 
0 measured from the upstream pole. All fluid equations are written nondimcnsionally with the radial 
coordinate scaled by a, and the velocity by U, the terminal velocity. Because the fluid motion in the 
surrounding liquid is axisymmctric and incompressible, the velocities in the r (u r ) and 0 (uo) directions can 
be specified in terms of a (nondimensional) stream function M ; (r,0) (sealed by a 2 U): 

I dy 1 dy 


u 


r sin# d6 


"« = 


/•sin# 3 


and the one component of the vorticity , w (in the azimuthal direction, sealed by U/a) is 

1 dy I 8 [ 1 dy 

rw(r,6) = ^+ — — 

sin# 3 % r 2 dO sin# 89 


(3) 


(4) 


The fluid motion is governed by the Navicr Stokes equations at order one Reynolds number, and we use the 
vorticity-strcam function formulation to develop solutions: 
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where Re is the Reynolds number (— pUa J\i). The boundary conditions on the stream function and vorticity 
are the matching to the free stream velocity at infinity (r— *») 


¥ -» 


— r^sinO 
2 


(6a,b) 


H — > ()j 

zero normal velocity, and the balance of the tangential stress with the Marangoni stress at the bubble 
interface r=l: 

i// = 0 


2 d¥ Ma ft (7a,b) 

Mi = — - H 

sin 9 ft 1 - f d6 


where T is the surface concentration (nondimensionalized by !"«,) and Ma is the Marangoni number RTIf* 
/pU. The surface concentration is obtained from the solution of the convective diffusion equation, 

m«VC = — V 2 C (8) 

Pe 

where C(r,9) is the bulk concentration (nondimensionalized by the far field value Co) and Pc the Peclct 
number, and the surfactant boundary conditions at the bubble surface r=l : 


r = 


k£s 

i + *c M 


J d_ 

sin# dO 


(sin 5T« fl ) = — — 
e Pe dr 


(9a,b) 


where k=Co/b and x = ba/Too as in Sec I. The nondimensional concentration must match to 1 as r->x 

Wc use a finite difference method to solve the above equations. The infinite boundary is truncated 
to a value r =20, and the annular region from the bubble surface r=l to the outer boundary r=20 is mapped 
to a unit square (x,y) by the transformations x=lnr/lnr a) and y= 7c/0. The square is discretized with Ax=.01 
and Ay=.01. The discretized system is solved by using an ADI (Alternating Directions Implicit) scheme 16 
in which the field equations (4,5 and 8) and boundary condition (9b) arc made pseudo-unsteady in time, and 
equal part time steps (At=.005), implicit in either x or y, are taken. The convergence criteria in time is 
formulated in terms of the change in the maximum value of the vorticity and the stream function over At 
10 2 ; this change is required to be less than 10' 6 

III Simulation Results 

We begin by first examining the case without inertia (Rc=0). To illustrate remobilization, wc fix 
the values of the Marangoni and Peclct numbers (Ma=5, Pc=l) and % (=1), and wc obtain solutions for 
increasing bulk concentrations by vary ing k. In Fig. 1, wc plot the surface concentration profile for three 
different values of k (=.1,1 and 5) relative to the maximum packing concentration (Fig. la) and the 
equilibrium concentration T c> (Fig. lb) The convective pardoning of surfactant on the surface is evident for 
all values of k, as the surface concentration is higher at the downstream pole (0=7c) then the upstream pole 
(0=0). We note that as k increases. Fig la indicates clearly that the total amount of adsorbed surfactant 
increases. More importantly, as is evident in Fig. lb, the aurfacc concentration becomes more uniform as k 
increases: The ratio of bulk diffusion to convection (x(l+k)/Pc increases from 11 to 6 as k v aries from . 1 
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to 5, and diffusion therefore begins to outscalc convection (see criteria 2b) to maintain a uniform surface 
concentration The remobilization of the interfacial mobility is evident in Figs. 2 which plots the maximum 
interfacial velocity as a function of k over a much wider range of k (10 2 <k<10 2 ) with the same values of 
Ma, Pc and x as in Fig 1 . Note that for a clean interface, the surface velocity is equal to (l/2)sin0, so the 
maximum velocity at the equator (0 =tc/ 2) is 1/2. As k increases from 0 to 1, the interfacial velocity 
decreases from the clean value of 1/2. This is the retardation that as we noted in the Introduction has been 
well documented in the literature. However, as k increases further and the surface concentration becomes 
uniform, the interfacial mobility is restored and the maximum surface velocity tends to .5 . 

As the bulk concentration increases and the surface concentration becomes more uniform, the 
sublayer concentration (in equilibrium with the surface) also becomes more uniform. In Fig. 2a, we plot the 
sublayer concentration for k= 1,1 and 5 (for Ma=5,Pc=l and x = l), ami the tendency towards a more 
uniform concentration with increasing k is evident. With the sublayer concentration approaching 1, the 
concentration in the bulk also approaches one as the diffusion driving force disappears. In Fig 3, contours 
of constant bulk concentration arc shown for k= 1 (Fig. 3a), k=l (Fig. 3b) and k= 1 0 (Fig. 3c), and it is clear 
that as the interface remobilizes, the bulk concentration becomes uniform. 

The effect of increasing the bulk concentration on the terminal velocity is examined in Fig. 4 which 
plots the drag on the nondimcnsional drag on the bubble (nondimcnsionalized by 7C|iaU) as a function of k, 
for Ma=5 and x=l and a few values of Pe (. 1,1 and 10). In these nondimcnsional units, the drag on a clean 
bubble is 4 and the drag on a completely immobile surface (the Stokes drag) is equal to 6. For fixed Pc, we 
note that as k increases, the drag at first increases (corresponding to the decrease in interfacial mobility 
observed in Fig. 2b), but then decreases as the interface becomes remobilized. For fixed k, as Pc increases 
the drag increases. This elevation in drag, which begins to approach the Stokes value of 6 for Pe=10, can 
be directly attributed to the larger convective partitioning of surfactant on the surface relative to bulk 
diffusion; as Pe increases, x( l+k)/Pe decreases. Importantly, the larger drags at high Pcclct numbers can 
also be reduced to clean surface drags if the bulk concentration is taken large enough so as to increase the 
ratio x(l+k)/Pc. 

The above simulations have been for Rc=0; some preliminary results for a finite value of Re 
(Rc=5) arc given in Fig. 5 in which the drag is plotted as a function of k for Pc=l and Ma-5. Again we 
note that as k increases, the drag can be reduced and the terminal velocity' increased. 
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Fig. 1 Surface concentration distribution as a function of k for Rc-0, Ma-5, x~l and Pc-1. 



Fig. 2. The maximum interfacial velocity as a function of k and the sublayer concentration for Re=0, 
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Fig. 4. Terminal velocity as a function of k for Re=0 



Fig. 5 Terminal velocity as a function of k for Rc=5 






